Scaling and universality in the 2D Ising model with a magnetic field 
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The scaling function of the 2D Ising model in a magnetic field on the square and triangular 
lattices is obtained numerically via Baxter's variational corner transfer matrix approach. The use of 
the Aharony-Fisher non-linear scaling variables allowed us to perform calculations sufficiently away 
from the critical point to obtain very high precision data, which convincingly confirm all predictions 
of the scaling and universality hypotheses. The results are in excellent agreement with the field 
theory calculations of Fonseca and Zamolodchikov as well as with many previously known exact and 
numerical results for the 2D Ising model. This includes excellent agreement with the classic analytic 
results for the magnetic susceptibility by Barouch, McCoy, Tracy and Wu, recently enhanced by 
Orrick, Nickel, Guttmann and Perk. 
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The principles of scaling and universality (see, e.g., 
play important roles in the theory of phase transition and 
critical phenomena. The scaling assumption asserts that 
observable quantities exhibit power law singularities in 
the variable AT = T — T c in the vicinity of the critical 
temperature T c , with coefficients being functions of cer- 
tain dimensionless combinations of available parameters, 
e.g., the magnetic field H and AT. The universality hy- 
pothesis states that the leading singular part of the free 
energy is a universal scaling function which is same for all 
systems in a given "universality class" . In two dimensions 
classes of universal critical behaviour are well understood 
— they are classified by conformal field theory (CFT) Q. 
The latter provides exact solutions for 2D systems at the 
critical point and, in particular, gives exact values of the 
scaling dimensions. However, the calculation of scaling 
functions describing off-critical behaviour is a hard prob- 
lem which does not have any exact analytical solutions. 
From the field theory point of view this requires solving a 
CFT perturbed by at least two relevant operators. Direct 
numerical calculations and simulations in lattice models 
also face serious difficulties due to multiple long range 
correlations and the resulting poor convergence near the 
critical point. 

It appears that despite numerous analytical and nu- 
merical results (cited below), the full picture of scal- 
ing and universality has never been convincingly demon- 
strated through numerical calculations in lattice mod- 
els. Our aim is to do this. Here we consider the planar 
nearest-neighbour Ising model on the regular square and 
triangular lattices, which has already played a prominent 
role in the development of the theory of phase transition 
and critical phenomena 2-8]. Its partition function reads 

z = j2^p{pY, a ^ + H J2^} ^ °-i = ±l . (!) 

a {ij) i 

where the first sum in the exponent is taken over all 
edges, the second over all sites and the outer sum over 



all spin configurations {a} of the lattice. The constants 
H and /3 denote the (suitably normalized) magnetic field 
and inverse temperature. The specific free energy, mag- 
netization and magnetic susceptibility are defined as 
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F = — lim 
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N"°~' " dH' dH 2 
where N is the number of lattice sites. The model ex- 
hibits a second order phase transition at H = and 
/3 = j3 c , where 

/3W=ilog(l + V2), /3« = ±log3, (3) 

for the square Q and triangular 0] lattices, respectively. 

The scaling and universality hypotheses predict that 
the leading singular part, F sing (A(3, H), of the free en- 
ergy in the vicinity of the critical point, A/3 = /3 — (3 C ~ 0, 
H ~ 0, can be expressed through a universal function 
T(m, h), 



F, 



smg (A[3, H) = F(m(A{3, H),h(A/3, H)), (4) 



where A/3 and H enter the rhs only through non-linear 
scaling variables , 

0(A/3) + O ((A/3) 3 ) + 0{H 2 ) + ..., 



m = m(A/3, H 
h = h(A(3,H)=0(H)+HO(Ap)+0{H 3 ) 



(5) 



which are analytic functions of A/3 and H . The coeffi- 
cients in these expansions depend on the details of the 
microscopic interaction (for instance they are different 
for the square and triangular lattices), but the function 
J-(m, h) is the same for all models in the 2D Ising model 
universality class. It can be written as 
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where $(17) is a universal scaling function of a single vari- 
able rj (the scaling parameter), normalized such that 
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The function T(m, h) has a concise interpretation in 
terms of 2D Euclidean quantum field theory. Namely, 
it coincides with the vacuum energy density of the "Ising 
Field Theory" (IFT) [llj. The latter is defined as a model 
of perturbed conformal field theory with the action 



A 



IFT 



■4(c=l/2) + 
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e(x) d 2 x + h / a(x)d 2 x, (8) 



where -4( c= i/2) stands for the action of the c = 1/2 CFT 
of free massless Majorana fermions, <j(x) and e(x) are pri- 
mary fields of conformal dimensions 1/16 and 1/2. The 
parameters m and h have the mass dimensions 1 and 
15/8, respectively, and the scaling parameter 77 in © is 
dimensionless. 

The scaling function © is of much interest as it con- 
trols all thermodynamic properties of the Ising model in 
the critical domain. Although there are many exact re- 
sults (obtained through exact solutions of © at h = 
and all r |, fl [HQ, and at r = and all h H HHHI; 
these data are collected in 22J) as well as much numeri- 
cal data [23- 28] about this function, its complete analytic 
characterization is still lacking. 

Recently ll| the function ^ , particularly its analytic 
properties, have been thoroughly studied in the frame- 
work of the IFT (T8T) . The authors of 111 made extensive 



The authors of 
numerical calculations of the scaling function $(77) for 
real and complex values of 77. 

One of the motivations of our work was to confirm 
and extend the field theory results of [ll[ through ab 
initio calculations, directly from the original lattice for- 
mulation (TTJ) of the Ising model. Here we give a brief 
summary of our results for the triangular and square lat- 
tices (the latter were previously reported in [29]). We 
used Baxter's variational approach based on the corner 
transfer matrix method [30I l3lj . The main advantage 
of this approach over other numerical schemes (e.g., the 
row-to-row transfer matrix method) is that it is formu- 
lated directly in the limit of an infinite lattice. Its accu- 
racy depends on the magnitude of truncated eigenvalues 
of the corner transfer matrix (which is at our control), 
rather than the size of the lattice. The original Bax- 



ter approach [31[ was enhanced by an improved iteration 
scheme 32|, known as the corner transfer matrix renor- 
malization group (CTMRG). The use of the non-linear 
scaling variables ([5]) allowed us to perform calculations 
sufficiently away from the critical point with a reliable 
convergence of the algorithm. The results for the scaling 
function $(r]) are shown in Fig. [TJ In total we have cal- 
culated about 10,000 data points for different values of 
the temperature and magnetic field on the square lattice 
and about 5,000 for the triangular one. As seen from 
the picture all points collapse on a smooth curve, shown 
by the solid line (as expected, the curve is the same for 
the square and triangular lattices). The spread of the 
points at any fixed value 77 does not exceed 10~ 6 relative 
accuracy. This gives a convincing demonstration of the 
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FIG. 1: The scaling function $(77) in the three regions, sepa- 
rated by the dashed lines at n ~ ±2.3, can be approximated 
with high precision (up to 10 -6 ) by the series (|10|) . (|lip and 
(£9j) with coefficients given in Table 1. 



scaling and universality in the 2D Ising model. Further- 
more, our numerical results for $(77) remarkably confirm 
the field theory calculations [llj], to within all six signif- 
icant digits presented therein. 

For further reference we write down expansions of the 
function $(77) for large values of 77 on the real line 



lSfc/8 



77 +00 , (9) 
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$ hig h(v) = ?7 2 E G ^r 30fe/8 , (10) 
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and for small values of 77, 



®(V) = -5- log?? 2 
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Some of the above expansion coefficients are known ex- 
actly. The coefficient Gi has a simple explicit expression 
5j; the coefficients G2 and G2 have integral expressions 
13 . [l3l ] involving solutions of the Painleve III equation. 
They were numerically evaluated to very high precision 
(50 digits) in [33|. The coefficients <&o,®i were analyt- 
ically calculated in Q and [34|, respectively. The nu- 
merical value of $1 (which requires certain quadratures) 
was found in [29j]. The above values are quoted in the 
last column of Table 1 . 

In what follows we exclude the temperature variable /3 
in favour of a new variable 



(1 -sinh 2 2/3)/(2sinh2/3), 



(sq. lat. 



\ {e-P -e? sinh 2/3) /(sinh 2/3) 1 / 2 , (tr. lat.) 

(12) 

which is vanishing for (3 — j3 c and positive for /3 < f3 c 
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(above the critical temperature). Another useful variable 

f 16e 8 ' 3 /(e 4 ' 3 - l) 4 , (sq. lat.) 

fc 2 = < (13) 
\ 16e 4 ^/((e 4 ^-l) 3 (e 40 + 3)), (tr. lat.) 

The lattice free energy for r, H — >• 0, 

F(t,H) =F aing (T,H)+F reg (T,H)+F aub (T,H), (14) 

contains leading universal part Q, regular terms 
Freg{T, i?), which are analytic in r and H, and sublead- 
ing singular terms F su b(r, H), which are non-analytic, 
but less singular than the first term in (| 14|) . Therefore, to 
extract the universal scaling function from the lattice cal- 
culations one should be able to isolate and subtract these 
extra terms. Moreover, one needs to know the explicit 
form of the non- linear scaling variables ([5]). In princi- 
ple, all this information can be determined entirely from 
numerical calculations (provided one assumes the values 
of exponents of the subleading terms, predicted by the 
analysis [H, 35 1 of the CFT irrelevant operators, con- 
tributing to the free energy HJ])). Much more accurate 
results can be obtained if the numerical work is com- 
bined with known exact results. Namely, the zero-field 
free energy reads 0, 0] 
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F( s ~) (t, 0) = — g log(4 sinh 2/3) — g^j- // dfadfa 

o 



log(2vT 



COS <p\ — COS <p2) 
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(15) 



FM(r,0) = -± log(4 sinh 2/3)-^ J J dfadfa 

o 

l0g(3 + T 2 - COS (j)i - COS 02 - cos(0i + 4> 2 )), 



where the superscipts (s) and (t) stand for the square 
and triangular lattices, respectively. Write the non-linear 
variables ([5]) in the form, 



m(r, H) 



h(r, H) 



-C t t a(T)+H 2 b(T) + 0(H 4 ), 



C h H 



c(r) + H 2 d(r) + 0(H 4 ) , (16) 



where a(0) = c(0) = 1, h(r,H) = —h(T,—H) and write 
the regular part in (fT4"]) as, 



F reg (r, H) = A(t) + H 2 B{t) + 0(H 4 ) 



(17) 



As shown in [33j , the most singular subleading term, con- 
tributing to (jT4")) is of the order of t 9 / 4 H 2 ~ m 6 for the 
square lattice and r 13 / 4 H 2 ~ m 8 for the triangular lat- 
tice. 

Rewriting (|15j) in the form (|14[) plus regular terms, one 
obtains 



and 



a W (T) = l_^ r 2 + ^ T 4 + 0(r 6 ); 



(18) 



(19) 



The contribution to the regular part reads 

_^ T 3 + 5(l^) r 4 + 0(T5); 



^)(r) = -ACl 2 (f)-ilog; 

l_\ T 2 7_ T 3 , f 4+9 log 12 1_\ 4 



2+3 log 12 

87TV3 



(20) 
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+ ( 4+9 log 12 _ J-) T 4 + 0(t 5 ) 



Next, with the definition (fT3|) the zero- field spontaneous 
magnetization has the same expression for both lattices 

M(r,0) = (1 - fc 2 ) 1 / 8 , r<0. (21) 

Combining this with (|2J|, ©, © and (JUJ one obtains 



2 3 / 16 /G!, Cf = -2 B / 16 3- 3 / 32 /Gi, (22) 



and 



(sU , , r 15t 2 9r 3 4333r 4 ^. 5 . 
c()(r) = 1+ 4 + W-512-^8304 +0(T) ' 



c (t) (r) = l + - + 



5r 2 



727r 4 



6 96 576 165888 



■0(r 5 



(23) 

Finally, consider the zero-field susceptibility. The second 
field derivative of (|14[) at H = gives 



X(r) 



2GG 2 c(r) 2 d 2 F stifc 



2B(t) 



(V2|r| a (r)) 7 /4 dff 2 
r a(r) b(r 



H=0 
2. 



(24) 



V2 7T 



l + log(2r 2 a(r))), 



where G — G 2 for r > and G = G 2 for r < 0. No 
simple closed form expression for the zero-field suscepti- 
bility x(t) is known. However, the authors of [33| ob- 
tained remarkable asymptotic expansions of x( T ) f° r the 
square lattice for small r to within 0(r 14 ) terms with 
high-precision numerical coefficients. Using their results 
in (|2~4"|) . one obtains 

£ (s) (t) = 0.0520666225469 + 0.0769120341893 r 



0.0360200462309 r 2 + 0(r 3 ), 



and 



&( s )(t) = /ijf 3 (l + I + 0(r 2 )) , = 0.071868670814 

(26) 

No similar expansion for r ^ is available for the trian- 
gular lattice. We used our data for r = to estimate 

B (t) (r) = 0.0247805582(2) + O(r), ^ = -0.010475(1) 

(27) 

and the coefficient d(r) = + O(t) in (|16l) 

e W = -0.00728(30), = +0.00129(1), (28) 
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which is in agreement with 



» 



-0.00727(15) from 



The above expressions were used to analyze our exten- 
sive numerical data and extract the necessary informa- 
tion to obtain the universal scaling function. The results 
are summarized in Table 1. For convenience of compari- 
son we quoted the field theory results from llj. Earlier 
exact and numerical results for the same quantities are 
also quoted (whenever available). 

To conclude, we have implemented Baxter's variational 
corner transfer matrix approach to obtain the universal 
scaling function for the Ising model in a magnetic field on 
the square and triangular lattice, as shown in Fig. 1 and 
Table 1. The numerical data is seen to be in remarkable 
agreement with the field theory results obtained by Fon- 
seca and Zamolodchikov llj. We also report a remark- 



able agreement (11 to 14 digits) between our numerical 
values for Gi , G2 and G2 and the classic exact results of 
Barouch, McCoy, Tracy and Wu [!, [H, [HI and a sim- 
ilar agreement between the values <E>o and $1 and the 
exact predictions 15|, [3J] of Zamolodchikov's integrable 
Eg field theory [§]. Interestingly, this Eg symmetry has 
now been observed in experiments on the transverse Ising 
chain [36| . 
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TABLE I: Numerical values of the coefficients G n , G„, $ 
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